Coarsening in a Driven Ising Chain with Conserved Dynamics 
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We study the low-temperature coarsening of an Ising chain subject to spin-exchange dynamics and 
a small driving force. This dynamical system reduces to a domain diffusion process, in which 
entire domains undergo nearest-neighbor hopping, except for the shortest domains - dimers - which 
undergo long-range hopping. This system is characterized by two independent le ngth scales: the 
average domain length L(t) ~ i 1//2 and the average dimer hopping distance £(t) ~ U Lit) ~ t 1 ^ 4 . As 
a consequence of these two scales, the density Ck{t) of domains of length k does not obey scaling. 
This breakdown of scaling also leads to the density of short domains decaying as t~ 5,/4 , instead of 
the t~ 3 ^ 2 decay that would arise from pure domain diffusion. 
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I. INTRODUCTION 

The approach to equilibrium in isotropic systems which 
are quenched from a high-temperature homogeneous 
phase to a low-temperature two-phase region is now rel- 
atively well understood The basic feature of such 
systems is that they typically organize into a coarsen- 
ing mosaic of single-phase regions, with a characteristic 
length scale that grows as a power law in time. For driven 
systems, on the other hand, considerably less progress 
has been made in understanding the coarsening dynam- 
ics, although the stationary properties have been thor- 
oughly investigated ||. In the presence of driving, the 
physically relevant coarsening mechanisms are those with 
conserved order-parameter dynamics. This would be ap- 
propriate, for example, for treating the phase separation 
of binary liquids or binary alloys under the influence of 
gravity [Q. 

In this spirit, Cornell and Bray || recently studied 
the coarsening dynamics of an Ising chain which is en- 
dowed with conserved spin-exchange Kawasaki dynamics 
and which is also subject to a driving field which favors 
transport of up spins to the right and down spins to the 
left. They argue that in the limit of low temperature 
and weak field, the spin dynamics is equivalent to domain 
diffusion (DD) B. In this DD process, up domains hop 
rigidly by one lattice spacing to the right, and down do- 
mains hop by one lattice spacing to the left. Due to this 
nearest-neighbor hopping, small domains are progres- 
sively "squeezed out" and the adjacent neighboring do- 
mains coalesce. This random walk mechanism leads both 
to a reduction in the number of domains and an increase 
in their average length, L(t) ~ t 1 / 2 . Numerical evidence 
was also presented that the density of domains of length 
k obeys the scaling form Cfc(t) ~ (k/L 3 ) exp(— k 2 /L 2 ) 
||. This further implies a i~ 3 / 2 asymptotic decay for 
the density of domains of fixed length. 

The goal of this paper is to show that there is a subtle 
but crucial difference between the dynamics of individ- 
ual spins in the driven Ising chain with spin conserving 
dynamics and the DD process. The fundamental point 
is that for the shortest domains - dimers - the spin-level 



dynamics results in long-range dimer hopping, with the 
average jump length growing as \[L. In contrast, for 
the DD process, dimers necessarily jump to the next do- 
main boundary, that is, they jump a distance of order 
L. Therefore dimers disappear more slowly than in the 
DD process, and the overall coarsening rate is dimer con- 
trolled (DC). This DC system is characterized by two in- 
dependent length scales: (a) the average domain length 
L(t), which is still proportional to t 1 / 2 , as in the DD pro- 
cess, and (b) the average dimer hopping distance, which 
is proportional to t 1 ^ 4 . As a result of the two length 
scales, the density of domains of fixed length asymptot- 
ically decays as i -5 / 4 , instead of the t~ z / 2 decay of the 
DD process. Correspondingly, the domain length dis- 
tribution in the DC process does not obey conventional 
scaling in the small- length limit, although the domain 
length distributions for both the DD and DC processes 
are visually similar. 

In Sec. II, we define the spin dynamics precisely and 
describe the correspondences to the DD and DC pro- 
cesses. We also discuss their essential differences and 
show how the latter process provides a faithful descrip- 
tion of the microscopic spin dynamics. We then present 
heuristic arguments which suggest new kinetic behavior 
for the DC process. Simulation results which support 
our basic arguments are presented in Sec. III. In Sec. IV, 
we outline a perturbative approach, based on a matched 
asymptotic expansion, which accounts for the observed 
breakdown of scaling in the domain length distribution 
for vanishingly small minority fraction. Sec. V contains 
both a summary and a brief discussion of open issues. 



II. GEOMETRICAL PICTURE OF THE 
DYNAMICS 

The microscopic system is a chain of Ising spins with 
nearest-neighbor ferromagnetic interaction J. The chain 
is subject to spin-exchange dynamics, where the only pos- 
sible re-arrangement process is the exchange of two anti- 
parallel nearest-neighbor spins. Thus the magnetization 
is manifestly conserved (we use a magnetic terminology, 



I 



+- A = 


4 J -E 


(i), 


-+ A = 


4J + E 


(ii), 


++ A = 


-E 


(iii), 


+- A = 


-E 


(iv). 



although a system with conserved dynamics naturally ap- 
plies to an alloy) . The exchange occurs at a rate propor- 
tional to e~ A / T , where A is the energy difference between 
the initial and final states, and T is the temperature (with 
Boltzmann constant set to unity) . There is also a driving 
field E which favors motion of up spins to the right and 
down spins to the left. The spin-flip events are: 
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The first two processes occur on domain boundaries, 
while the last two account for the motion of a single spin 
which is inside a domain of the opposite sign. The "for- 
ward" processes involve the energy change A, while the 
"backward" processes have energy change — A. 

Interesting dynamics arises in the limit of low (but 
non-zero) temperature and weak driving field, that is, 
< T <C E <C J. To appreciate the nature of the dynam- 
ics for this parameter range, notice that in one dimension 
(ID) the order-disorder transition occurs at T c — 0. At 
T = 0, the spin-exchange dynamics traps the system in 
a metastable state which consists of domains of lengths 
> 2 HfjJ . To avoid this "freezing" , the temperature must 
be non-zero. At low but non-zero temperature, the sys- 
tem will coarsen as long as the mean domain length is 
smaller than the correlation length £ ~ e J l T 1. 

The limit where the driving field satisfies T <C E -C J 
leads to an approximate equivalence with the DD process 
||. To understand this correspondence, consider the sit- 
uation after domains have coarsened to relatively large 
length. By process (i — s -), an up spin may detach from 
the right edge of an up domain with rate e~\^ 1 ~ E )l T , or 
equivalently, a down spin may detach from the left edge 
of a down domain. Similarly, an up spin may also de- 
tach from the left edge of an up domain (or a down spin 
may detach from the right edge of down domain) by step 
(ii — i ). However, this process occurs at a rate which is 
a factor e~ 2E l T smaller than (i — x ). Moreover, even if 
step (ii — occurs, the detached spin quickly recombines 
with the same domain by the reverse process (ii t— ), since 
the motion of the detached spin away from the domain 
is energetically unfavorable. 

Once (i — v ) occurs, the system evolves further either 
by (iv — which corresponds to the up spin moving to 
the right and eventually joining the next up domain, or 
by (iii — *•), where a down spin moves to the left and joins 
the next down domain. The former process is illustrated 
in Fig. 1. As a result of these processes, an up domain 
hops rigidly by one lattice spacing to the right or a down 
domain hops one spacing to the left. If wc measure the 
time in units of e^ J ~ B ^' T , then both these hopping steps 
occur at unit rate, independent of the domain length. 
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FIG. 1. Illustration of the detachment of an up spin from 
an up domain and its merging with the neighboring up do- 
main to the right. In the upper part of the figure, each line 
represents the state of the system after a single spin-exchange 
event. This evolution is equivalent to rigid-body domain hop- 
ping, with a down domain hopping one lattice site to the left, 
as indicated in the lower portion. 

An essential feature of the low-temperature weak-field 
limit T -C E -C J is that all other processes are asymp- 
totically negligible in the intermediate-time range where 
coarsening is occurring. Thus in this time range the sys- 
tem consists of the contiguous array of alternating up 
and down domains, and the dynamics proceeds by tak- 
ing an up (down) domain and moving it one lattice site to 
the right (left). Whenever a domain shrinks to zero size, 
its two adjacent neighbors coalesce. This description is 
the basis of the correspondence to the DD model. It is 
also worth noting that in the absence of a driving field, 
the dynamics again reduces to a DD process, but with 
a length-dependent hopping rate that is proportional to 
the inverse domain length . 
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FIG. 2. Time evolution of an up dimer in a sea of down 
spins. When the dimer dissociates, the isolated spins inde- 
pendently hop to the right. Shown is the space-time trajec- 
tory of the dimer for the case where the dimer recombines 
and becomes stationary again before the next domain wall is 
reached. 

A crucial feature of the mapping between the spin and 
the domain dynamics which is not apparent from the 
above description is the evolution of dimers (see Fig. 2). 
Consider an up dimer. If the rightmost spin of the dimer 
detaches, the dimer is converted into two isolated up 
spins in a sea of down spins. According to the spin dy- 
namics, each isolated spin independently and freely hops 
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to the right. Consequently their separation undergoes a 
simple random walk. The motion of spins anti-parallel to 
the field can be neglected, since this motion is inhibited 
by a factor e~ 2E l T . The hopping of this pair of separated 
spins terminates in one of two ways: (i) The rightmost 
up spin reaches the next domain boundary and subse- 
quently the other up spin hits this same boundary. This 
corresponds to the coalescence of the two adjacent down 
domains and is part of the DD picture, (ii) The dimer 
recombines before the next domain boundary is reached 
(Fig- 1). 

Dimer recombination is the crucial new feature which 
was not included in the DD process. This recombination 
plays an essential asymptotic role because the average 
dimer jump distance I is much smaller than the average 
domain length L in the long-time limit. Consequently, 
recombination of a dimer is much more probable than 
domain coalescence. To verify this assertion let us esti- 
mate I. The dimer recombines if the separation between 
the two spins shrinks to zero before the rightmost spin 
reaches the next domain boundary. This is a classic first- 
passage process, and the probability that this separation 
first equals zero at 21 steps is given by || 

va ) = J_ («-2)' (1) 

l() 2«-i(I-l)!i! [) 

For large I, this expression simplifies to V(l) ~ 
7T -1 / 2 ^ 3 / 2 . The average dimer jump distance £ may now 
be estimated as I = X)z<l ~ V^- Thus asymptot- 

ically i < L. 

We now use this picture to estimate the overall time 
dependence of small-length domains. The crucial feature 
is that a domain can disappear only if a dimer first disso- 
ciates and then does not recombine before its constituent 
spins reach the next domain boundary. From the analogy 
with the first-passage process, the probability R(t) that 
the dissociated dimer does not recombine before the next 
domain boundary is given by 

R-jrv^-jrr^-L- 1 / 2 . (2) 

l=L l=L 

Since the disappearance of a dimer leads to domain co- 
alescence, the total number of domains N(t) obeys the 
rate equation 

dN ur ° 2 ™ 

On the other hand, the dynamics of large domains 
should still be governed by the gain and loss of single 
spins at the boundary, as outlined in Fig. [|. Since these 
gain and loss processes occur at the same rate, the do- 
main length undergoes an isotropic random walk, so that 
L(t) should grow as t 1 ^ 2 . Correspondingly, the number of 
domains N(t) decays as t -1 / 2 , the inverse of the average 
domain length. Substituting these two expectations into 



Eq. (||), we immediately obtain 

C 2 {t)~t- b l\ (4) 

This slower decay of dimers is one of the primary fea- 
tures of the DC process. It should be compared with 
the prediction C%(f) ~ £~ 3 / 2 which arises from the DD 
process. This latter time dependence would be obtained 
from Eq. (||) if the rate of dimer disappearance were 
unity, rather than proportional to L~ 1//2 . 

This —5/4 exponent value has fundamental implica- 
tions for the density of domains of length k, Ck(t). Let 
us suppose that this density obeys the conventional scal- 
ing hypothesis 

Cu{t) ~ > (|) , (5) 

where the prefactor L~ 2 follows from the length normal- 
ization condition J2 kC k(t) = 1. If C 2 (t) ~ i~ 5/4 , then 
either F(z) ~ \fz as z — * 0, or Ck(t) does not obey scal- 
ing for small k. We shall present evidence from both 
simulations and an analytical approach that strongly fa- 
vors the latter alternative. 



III. SIMULATION RESULTS 

In our simulations, we first initialize an array of al- 
ternating up and down domains of random lengths. For 
minority phase density fj,, we choose the average length of 
minority domains to be L = 10, and L/ jj, for the majority 
phase. The time evolution involves the following steps: 
(i) Pick a domain at random, (ii) Move an up (down) do- 
main of length > 2 to the right (left), (iii) If the domain 
is a dimer, choose its jump distance I from the probabil- 
ity distribution V(l) given by Eq. (|l|). If I exceeds the 
length of the neighboring domain, remove the dimer and 
merge the surrounding domains, (iv) Update the time 
by 1/ (number of domains). Simulations were performed 
on a chain of 4 x 10 6 domains for times up to 5 x 10 5 , 
and averaged over 16 samples. This is of the same order 
of data as the simulations of Cornell and Bray . 



A. Average time-dependent properties 

In Fig. ||, we plot N(t) for various minority fractions 
\x. Also shown are the corresponding results for the DD 
process. The linearity of the data suggests power-law be- 
havior, and visually the asymptotic slopes are very close 
to the expected value of —1/2. The essential difference 
between DC and DD models is manifested by the be- 
havior of the dimer density. Fig. ^ shows that the dimer 
density indeed decays more slowly than in the DD model. 
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FIG. 3. Time dependence of the domain density for vari- 
ous minority spin fractions fi. Upper set of points - DC data, 
lower set - DD data. The DD data are divided by 2 to sepa- 
rate the two sets. As a guide to the eye, the solid line has a 
slope —1/2. 




FIG. 4. The density of dimers as a function of time. Upper 
points - DC, lower points - DD. As a guide to the eye, the 
solid lines have slopes —5/4 (upper) and —3/2 (lower). 

To highlight this difference, we plot the corresponding 
local exponents in Fig. ||. We define the local exponent 
at time t as the best-fit straight line to 10 successive 
data points (equally spaced on a logarithmic scale) up to 
time t in the double logarithmic plot of G^if) versus t. 
This definition significantly smooths statistical fluctua- 
tions while still revealing systematic trends in the data. 
As shown in Fig. ||, the local exponents of the DD and 



DC models are clearly different. For the DD model, these 
exponents are close to the expected value —3/2 and also 
exhibit weak systematic time dependence. Thus, the nat- 
ural conclusion is that the dimer density (as well as the 
density of domains of any fixed length) asymptotically 
decays as i~ 3 / 2 . 



For the DC model, the situation is more subtle. The 
local exponents initially are increasing with time, but this 
time dependence slows when the effective exponent value 
is close to the anticipated value —5/4. However, the sys- 
tematic ambiguities in the data make an extrapolation for 
the asymptotic value of the exponent uncertain. This un- 
certainty and the relatively small difference in the dimer 
exponent for the two models led us to consider C 2 {N) 
rather than C^if)- Indeed, since both N(t) and C^if) 
should be influenced by the same pre-asymptotic correc- 
tions, such corrections might cancel when C 2 is expressed 
as a function of N. From a scaling perspective, it is also 
natural to express dependences in terms of intrinsic vari- 
ables, rather than in terms of the extrinsic time variable. 
Thus in Fig. we plot the local exponents of C 2 (N) ver- 
sus N and the results are now relatively straightforward 
to interpret. Since the variation in the local exponent for 
the DC model is small over the entire range of N, the 
result C 2 (N) <~ N 5 ^ 2 is strongly suggested. Similarly, 
the DD data suggests that C2(N) ~ N 3 , as is antici- 
pated from N(t) ~ t- 1 / 2 and C 2 (t) ~ t~ 3/2 . Coupled 
with the basic result N(t) ~ t~ x / 2 which holds for both 
models, we conclude, now with considerable confidence, 
that C 2 (t) ~ i" 5/4 for the DC model. 
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FIG. 5. Local exponents for the time dependence of the 
number of dimers. Upper points - DC model, lower points - 
DD model. 
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At the scale shown, this distribution appears to exhibit 
data collapse. In fact, at the resolution of this figure, 
the domain length distributions for the DD and DC pro- 
cesses are virtually indistinguishable. However, for the 
scaling form to be compatible with Cz(t) ~ i _5/4 , the 
scaling function must vary as F(z) ~ z 1/2 as z — > 0, 
while the length distribution appears to be linear in k 
in the small- length limit. This linearity implies that the 
distribution cannot obey single-parameter scaling for all 
lengths. In fact, a closer examination of the small-fc tail 
(Fig. |^) reveals that there is that there is a small but 
systematic deviation from data collapse. The breakdown 
of scaling is manifested by the length distribution having 
a non-zero intercept with the k = axis, whose value is 
systematically decreasing with time. 
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FIG. 6. Local exponents for the dependence of the number 
of dimers on the total number of domains. Lower points - 
DC model, upper points - DD model. 



B. Domain length distribution 

The behavior of the domain length distribution is es- 
pecially interesting, as scaling breaks down in the small- 
length limit. This ultimately arises from the multi-step 
dissociation and recombination processes that govern the 
disappearance of dimers. We first test the conventional 
scaling hypothesis for the domain length distribution, 
namely, C k {t) ~ L~ 2 F{k/L), by plotting L 2 C k (t) ver- 
sus k/L in Fig. fj]. 
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FIG. 7. Scaling plot for the domain-length distribution for 
equal fractions of up and down spins. Notice that this distri- 
bution does not reach zero at k/L = (see Fig. 0). 



FIG. 8. Small-length tail of the domain length distribution 
for equal fractions of up and down spins, showing the system- 
atic departure from scaling and an intercept at k — which 
is systematically decreasing with time. 



IV. DOMAIN LENGTH DISTRIBUTION IN THE 
LIMIT n 

To better understand the nature of the domain length 
distribution, we focus on the limit where the fraction of 
minority spins /i is vanishingly small. This leads to con- 
siderable simplification in the domain dynamics. Gen- 
erally, the length of a domain can change by ±1 due 
to diffusion of neighboring domains, or the length can 
change arbitrarily by domain coalescence. In the limit 
fi — > 0, the diffusive "shrinkage" governs the disappear- 
ance of minority domains, while coalescence governs the 
disappearance majority domains. To verify this, let us 
estimate the characteristic times for the disappearance 
of the majority and minority domains by shrinkage. Let 
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L_ (£+) denote the average length of minority (major- 
ity) domains. A majority domain can shrink to zero in 
a time t+ of order L%, while the shrinking and disap- 
pearance of a minority domain requires a time t— ~ L 2 _ . 
Thus 



t- 



17 



(6) 



Therefore in the minority limit, shrinkage of majority do- 
mains, or equivalently, coalescence of minority domains 
is negligible. Consequently the minority domains are ef- 
fectively non-interacting and they evolve only by the ad- 
dition or loss of single spins as a result of the hopping 
of majority domains. Therefore the density of minority 
domains Ck(t) obeys the discrete diffusion equation 



dCk 
dt 



= Cfe+i - 2C k + Ck-u k>2. 



(7) 



The density of dimers (k — 2) obeys a separate, but 
similar equation. For dimers, there is no gain term due to 
processes which involve monomers, and the loss of dimers 
due to their dissociation into two monomers and ultimate 
domain coalescence occurs at a rate R ~ L~ x l 2 (see Sec. 
II), since a dissociated dimer may recombine before the 
coalescence occurs. Therefore the master equation for 
C 2 (t) is 



dC 2 
dt 



Ca — Co 



C2_ 

yfL 



In the continuum limit, Eq. (]?]) is equivalent to 
dC k (t) d 2 C k {t) 



dt 



dk 2 



(8) 



(9) 



while Eq. (jg) provides the boundary condition. In this 
equation, the left-hand side scales as t~ 1 C2, while the 
last term on right-hand side scales as L~ X / 2 C2 ~ t~ 1 / 4 C2- 
Thus in the long-time limit the left-hand side is negligible 
and Eq. (0) becomes 



8C k (t) C k (t) 



dk 



0. 



(10) 



This radiation boundary condition |T(]] expresses the fact 
that a dimer does not necessarily disappear when it dis- 
sociates, but it may be reconstituted and then grow into 
a finite-size domain. 

By dimensional analysis, Eq. (^|) implies the existence 
of the usual diffusive length scale L — \/t. By similar rea- 
soning, Eq. (|l(]) suggests the existence of an additional 
length scale t — \[~L. The competition between these two 
scales determines the asymptotic behavior. We therefore 
separately consider the "inner" region of small domains 
k -C L and the "outer" region of large domains k £, 
and then match these limiting solutions in the overlap 
region I < k < L 111. 



In the inner region k <C L, the diffusion equation sim- 
plifies to ^p- = 0, whose solution is (7fc(t) = A(t)+B(t)k. 
Employing the boundary condition Eq. (BO) we obtain 



C k (t) inner = A(t) 1 + -t= 



k 

7l 



(ii) 



In the outer region k 3> £, the system is governed by the 
original diffusion equation (||), while Eq. ( [To| ) reduces to 
the absorbing boundary condition. The solution in this 
region thus becomes 



C k {tf 



t 3/2 



exp 



At 



(12) 



The inner and outer solutions should match in the over- 
lapping region < C I « L. This determines the ampli- 
tude A(t) in Eq. (|ll|) to be proportional to t~ 5 / 4 . The 
inner solution now becomes 



/(\inncr T i ^ 



(13) 



with 7 a constant. 

These two limiting forms for Ck(t) match smoothly in 
the overlap region £ <C k <C L. This further suggests 
that the domain length distribution for the entire length 
range can be accounted for by the composite form 



ex P _ TT 



(14) 



To determine the validity of this hypotheses, we test for 
the existence of the kjt?l 2 correction term, since the lead- 
ing i~ 5 / 4 time dependence has already been established. 
For this purpose, consider C%{£) — G^t) versus t. This 
difference eliminates the leading i~ 5 / 4 behavior and thus 
isolates the kjt?l 2 correction term (Fig. ^). As seen in 
the figure, the data for C3 — C2 is consistent with a t^l 2 
time dependence. This test also supports the correctness 
of the composite form of Eq. ([lj]) for the domain length 
distribution. Finally, our numerical data suggests that 
this same dependence holds for all values of fi. 
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assumption of statistical independence, it is possible to 
solve the rate equations for Ck(t). This solution repro- 
duces the correct dynamical exponent, as well as the lin- 
ear small-length tail for the domain length distribution 
]l2| . This approach further predicts an exponential de- 
cay for the the large-length limit of the domain length 
distribution, for any non-zero fraction of minority spins. 
However, our numerical simulations at zero magnetiza- 
tion suggest that this large-length tail has the leading 
behavior exp(— (k/L) u ), with v greater than 1 and less 
than 2. This puzzling feature, neither an exponential nor 
a Gaussian decay for the large-length tail of the distribu- 
tion, deserves more careful attention. 

We gratefully acknowledge partial support from NSF 
grant DMR9632059 and ARO grant DAAH04-96- 1-0114. 



FIG. 9. Plot of Cz(t) - C2(t) versus time. As a guide to 
the eye, a straight line of slope —3/2 is shown. 

V. SUMMARY AND DISCUSSION 

We investigated the low-temperature coarsening of an 
Ising chain subject to spin-exchange dynamics and a 
weak driving force. The spin dynamics was reduced to 
a dimer-controlled domain diffusion process. From this 
picture, we established the existence of two growing char- 
acteristic length scales; one is the fundamental diffusion 
length i 1 / 2 , which provides the average domain size, and 
the other, I ~ t 1 ^ 4 , is the average dimer hopping dis- 
tance. The competition between these two scales leads to 
a breakdown of conventional scaling in the small-length 
tail of the domain distribution. As a consequence, the 
density of fixed- length domains decays as i~ 5 / 4 as t — ► oo. 
A key step to verify this latter result was to study the 
dependence of Ck on N, rather than the dependence on 
t. 

For the one-dimensional system, several basic unre- 
solved issues remain. Thus far, an analytical solution 
for the length distribution of minority domains only has 
been obtained in the extreme minority limit. It would 
be worthwhile to study analytically the case of arbitrary 
minority fraction. One approach is to treat domains as 
statistically independent. Such an approximation is ex- 
act in the extreme minority limit and also works well 
for the coarsening of the undriven Ising chain for both 
spin-flip and spin-exchange dynamics [p^|-p^| . Under the 
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